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Abstract 

■ We present the formalism for the covariant treatment of gravitational radiation in a 
magnetized environment and discuss the implications of the field for gravity waves in the 
cosmological context. Our geometrical approach brings to the fore the tension properties of 
the magnetic force lines and reveals their intricate interconnection to the spatial geometry of 
a magnetised spacetime. We show how the generic anisotropy of the field can act as a source 
of gravitational wave perturbations and how, depending on the spatial curvature distortion, 
the magnetic tension can boost or suppress waves passing through a magnetized region. 

(N 
O . 

1 Introduction 



(N 
> 
in 



X 



Magnetic fields appear everywhere in the universe. Planets, stars and galaxies carry fields that 
are large and extensive. The Milky Way has a magnetic field, coherent along the plane of 



(2jT), the galaxy, with strength of a few fiG. Magnetic fields also permeate the intracluster region 

of galaxy clusters, extending well beyond the core region of the cluster. In addition, reports of 
Faraday rotation in high redshift Lyman-a systems suggest that dynamically significant magnetic 
fields might have been present at high redshift condensations. In short, the more we look for 
magnetic fields in the universe the more ubiquitous we find them to be Gravitational 
waves, on the other hand, are as elusive as ever. An inevitable prediction of general relativity, 
gravitational waves are propagating ripples in the spacetime fabric triggered by changes in the 
matter distribution. Their extremely weak coupling to matter, however, means that detecting 
these ripples is a formidable task. Given the ubiquity of magnetic fields in the universe and 
their strong presence near some of the most promising candidates for detectable gravity wave 
signals, understanding the interaction between the two sources becomes especially interesting. 
Cosmology could provide the grounds for an exploratory first step, and recently there were two 
attempts in this direction [Q, || . The former employed the covariant formalism to investigate the 
evolution of gravitational waves in a magnetised cosmological environment. The latter adopted 
a metric based approach to study gravity wave production by stochastic large-scale magnetic 
fields]^] Here, following on the work of 0, we attempt to revisit the issue. Our geometrical 
approach brings to the fore the vector nature and the tension properties of magnetic fields and 
reveals the intricate interconnection between magnetic tension and spatial geometry. A direct 
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x The issue of magnetically induced tensor (i.e. gravitational wave) signatures in the Cosmic Microwave Back- 
ground spectrum has been raised and discussed in [^, [?]] and also in Q]. 
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implication of this coupling is that the magnetic contribution to spatial curvature deformations, 
along the field's own direction, is always zero. 

We begin with a brief outline of the covariant formalism and then apply it to the treatment 
of gravitational wave perturbations in a weakly magnetized, almost-FRW, cosmological environ- 
ment. Our background universe is a non-magnetized, spatially flat FRW cosmology filled with 
a single perfect fluid of very high conductivity. This model is perturbed by allowing for weak 
gravitational waves and a weak magnetic field and the impact of the field on the evolution of 
the waves is analysed. We introduce no a priori relation between the wave and the magnetic 
anisotropies, which enables us to look beyond the gravity-wave production issue. Throughout the 
analysis we employ covariantly defined tensors and use scalars that are invariantly constructed 
from these tensors. This means that we can account for the full spectrum of the directionally 
dependent magnetic effects, particularly those coming from the tension properties of the field 
lines. Confining to large scales we calculate the magnetic effects analytically both in the ra- 
diation and the dust dominated eras. We find that the presence of the field leaves the linear 
evolution of the waves unaffected, but modifies their magnitude. This implies that the overall 
magnetic impact depends entirely on the initial conditions. In the absence of gravity waves, we 
find that the anisotropic nature of the field leads to wave production. In general, however, the 
magnetic effect is particularly sensitive to the initial curvature deformation as measured along 
the direction of the field lines. This is where the subtle role of the magnetic tension becomes 
important. A negative curvature perturbation is found to boost the energy of the waves, whereas 
positive curvature leads to a damping effect. In either case the tension properties of the field 
tend to keep the spatial curvature deformation down to a minimum. 

The layout of the paper is as follows. In Sees. 2 and 3 we provide a brief description of the 
covariant treatment of cosmological perturbations and of cosmic magnetic fields respectively. 
Section 4 presents the necessary formalism for the study of gravitational waves within weakly 
magnetized almost-FRW cosmologies. The equations for the linear evolution of the wave's energy 
density are given in Sec. 5, and in Sec. 6 we look into the magnetic effects on gravitational 
radiation during the radiation and the dust dominated eras. 

We employ a metric with signature ( — h ++), the spacetime indices take the values 0, 1, 2, 3 
and use geometrised units with c = 1 = 8ttG. The Riemann and Ricci tensors are fixed by 
2V[ a V 6 ]ii c = R a bcdU d and R ab = R c acb respectively. 



2 Covariant description of cosmological perturbations 

The Ellis-Bruni approach to cosmological perturbations ||, follows earlier studies by Hawking || 
and Olson Jul, and is based on Ehlers' work on covariant hydrodynamics fill]- For a detailed 



and updated presentation of the covariant formalism the reader is referred to [12|. The essence of 
the method is to identify a set of covariantly defined variables that describe the inhomogeneity 
and the anisotropy of the universe in a transparent and gauge invariant manner. The non- 
linear equations governing the dynamics of these variables are then obtained from the full field 
equations. The transparency of the covariant variables means that the physical and geometrical 
content of their dynamical equations is simple to extract. The non- linear formulae can be 
linearised about a chosen background yielding a set of equations that describe deviations from 
inhomogeneity and isotropy in a straightforward way. 

In the covariant approach all physical and geometrical quantities are decomposed with respect 
to a fundamental timelike velocity field u a . The latter is determined by the motion of the matter 
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in the universe and introduces a unique time plus space (1+3) threading of the spacetime, as 
opposed to the 3+1 ADM slicing. Every observer has an instantaneous three-dimensional rest 
space, which is the tangent hypersurface orthogonal to u a . In a general cosmological model u a 
is chosen so that it reduces to the preferred velocity at the FRW limit, thus guaranteeing the 
gauge invariance of the formalism. The metric h a b = g a b + u a Ub of the observer's rest space 
(g a b is the spacetime metric) also defines the projected covariant derivative D a according to 
Da7b... c " = h a d hb e . . .h c j . . . VrfT e ../"\ When u a is irrotational, D a reduces to the covariant 
derivative in the hypersurfaces orthogonal to the observer's world line. It is also convenient to 
introduce the derivative Tb... c '" = u a V a Tb... c '" along the flow lines of u a . With these definitions, 
the covariant derivative of the fundamental velocity field decopmoses as 

VfeMa = U ab + U3 a b + l&Kb - U a U b . (1) 

The above defines the shear a a b = D/ 6 u tt \, the vorticity ui a b = Dp,u tt ], the (volume) expansion 
scalar G = D a u a and the 4-acceleration ii a = u 6 V&u a associated with the observer's motion.^ In- 
troducing the projected totally antisymmetric tensor e a b c = rj a b c dU d , where r\ a bcd is the spacetime 
alternating tensor, we define the vorticity vector u) a = e a b c 0J hc /2. The latter is also written as 
u> a = — curlu a /2 with curki a = e a bJ^ h u c . In addition we introduce the generalised curl of tensors 
by curlT a 6... c = e de ( a D d T b ... c ) e . The volume expansion defines an average length scale a, namely 
the scale factor of the universe, via a/ a = 0/3. The variables a a b, w a , and ii a characterise 
anisotropy in the local expansion and vanish identically in the FRW limit. The projected gradi- 
ents of scalars, vector and tensors describe local inhomogeneity in the observer's instantaneous 
rest space and also vanish in exact FRW spacetimes. 

An additional key decomposition is that of the matter stress-energy tensor T a b . Relative to 
a fundamental observer, 

T ab = [J-UaUb + ph ab + 2q^ a U b) + TTab , (2) 

where \i = T a bU a u b is the energy density, p = h a bT ab /3 is the isotropic pressure, q a = h a b Tb c u c 
is the energy flux and ir a b = T/ a y\ are the anisotropic stresses of the matter component. In 
the FRW limit q a and 7T a b vanish, as the energy-momentum tensor is always of the perfect-fluid 
form. 



3 The magnetic field 

3.1 Covariant description of cosmic magnetic fields 



The covariant description of electromagnetic fields was given in [14], and the Ellis-Bruni ap- 
proach was applied to magnetized cosmological perturbations in fTq , |l6| ]. Covariantly, the elec- 
tromagnetic Faraday tensor (Fab) decomposes into an electric (E a ) and a magnetic (H a ) field 
as 

F a b = U[ a Fb] + £ a bcFt c , (3) 



2 Round brackets indicate symmetrization and square ones antisymmetrisation. Angled brackets indicate the 
projected, symmetric, trace-free part (PSTF) of second rank tensors (i.e. S/ k m = [h^ a c ht) d — (l/S)h cd h a b]S a b) and 



the projected part of vectors (i.e. V< a > = h a b Vb) fl3[| . 
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implying that E a = F ab u a and H a = e abc F bc /2. Assuming that the universe is described by 
an infinitely conductive medium during most of its lifetime, we can ignore the presence of the 
electric field. Then Maxwell's equations reduce to one propagation equation 



H {a) =-l&H a + a ab H b + e abc H b io c , (4) 



and three constraints 



e abc u b H c + curliia = J {a) , (5) 
2oj a H a = p e , (6) 
V a H a = 0, (7) 

where = h a b J b is the projected current density and p e = —u a J a is the charge density [p|]. 
The above determine the evolution of the magnetic field completely. Relative to the fundamental 
observer, the stress-energy tensor of the field decomposes as 



T ab = \H 2 u a u b + \H 2 h ab + Ii ab , (8) 



where H 2 = H a H a and 



n ab = -H {a H b) = \H 2 h ab - H a H b . (9) 



Note that in the absence of electric fields the electromagnetic Poynting vector vanishes. Thus, 
the magnetic field behaves as a special imperfect fluid with energy density p m = H 2 /2, isotropic 
pressure p m = H 2 /6 and anisotropic pressure U ab . The latter reflects the vector nature of the 
field and carries the tension properties of the magnetic force lines. 

3.2 The magnetic tension 

Magnetic fields exert an isotropic pressure in all directions and carry a tension along their lines 



of force, with each flux tube behaving like an infinitely elastic rubber band [[L7], 18 1. The tension 
properties are encoded in the eigenvalues of H ab . Orthogonal to H a there are two positive 
eigenvalues equal to 1/3 each. Thus, the magnetic pressure perpendicular to the field lines is 
positive, reflecting their tendency to push each other apart. In the H a direction, however, the 
eigenvalue is -2/3 and the magnetic pressure is negative. The minus sign reflects the elasticity 
of the magnetic lines and their tendency to remain as "straight" as possible. It should be 
emphasized that the total magnetic pressure along the direction of the field lines is also negative 
and equals p mg = — /x mg = —H 2 /2. In other words, the false vacuum condition is satisfied along 
the magnetic lines of force. 

The magnetic tension has also non-trivial implications for the geometry of a magnetised 
3-dimensional space. Consider, for example, the non-linear Gauss-Codacci equation associated 
with a non-rotating magnetized spacetime. The Ricci curvature tensor of the spatial hypersur- 
faces decomposes irreducibly to [15| 

7Z ab = ^lZh ab + \li a b — \®&ab + V c (aV C b) + E ab , (10) 

where 

n = 2{p + \H 2 )-lQ 2 + 2o 2 , (11) 
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is the associated Ricci scalar and we have assumed a infinitely conductive perfect fluid for 
simplicity. Ignoring all sources but the magnetic field and then contracting twice along the field 
lines we obtain 

K = = \H 2 + ilWV = , (12) 

on using expression (0) for II a (,. Note that r] a = H a / V H 2 is the unit vector parallel to the 
magnetic force lines (i.e. rj a u a = and rj a rj a = 1). Hence, the tension properties of the field 
ensure that the spatial curvature along the magnetic direction is unaffected by the energy input 
of the field. This result demonstrates the generic tendency of the field lines to remain straight 
and it is indicative of what one might call a natural magnetic preference flat spatial geometry. It 
is also independent of the energy density of the field, namely of how close together or far apart 
the magnetic lines are. Thus, no matter how much energy one pumps into or removes from the 
field the magnetic contribution to TZ a i,r] a rj h remains zero. 

It should be emphasised that, despite their directional dependence, the tension properties of 
the field can affect average scalars such as the volume expansion of the universe (see |fT9|| - |21|| ), 
As it turns out, they can also affect the energy density of gravitational waves passing through a 
magnetised region. 



4 Covariant description of magnetized gravitational waves 
4.1 The background equations 

Consider an unperturbed non-magnetized FRW universe filled with a single barotropic fluid of 
infinite conductivity. When the spatial sections are flat, the background model is described by 
two propagation equations 

ft = -(i + w)e M , (13) 

6 = -l& 2 + ^(l + 3w), (14) 

and one constraint 

V = \®\ (15) 

where w = p/pb. Assuming that the fluid retains the barotropic equation of state (i.e. p = 
p(fi) always) and remains highly conductive, we perturb the background allowing for weak 
gravitational waves and a weak magnetic field. The infinite conductivity of the medium means 
that we can disregard any large-scale electric fields. The weakness of the magnetic field allows us 
to treat its energy density, its isotropic pressure and also the anisotropic magnetic stresses as first 
order perturbations. This automatically ensures that all three of them are gauge- independent 
variables. Finally, we assume that the magnetic field is coherent on all scales of interest. 



4.2 The linear equations 

The covariant description of gravitational waves, in the absence of magnetic fields, was originally 
considered by Hawking || , while more recent treatments can be found in 23 1 .[] Covariantly, 

3 Although the magnetic field can be treated as a viscous fluid, the approach discussed in |)j is not applicable 
here, since we are not introducing any phenomenological relation between the shear and the anisotropic stresses. 
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we monitor gravity waves via the electric (E ab = E^ ab \) and magnetic (H ab = H/ ab \) components 
of the Weyl (or conformal curvature) tensor C abcc [. The latter describes the locally free gravita- 
tional field, namely tidal forces and gravity waves. For a fundamental observer the Weyl tensor 
decomposes as 

Cabcd {dabqpdcdsr VabqpVcdsr) U^U E^ {j]abqp9cdsr 9abqpVcdsr) U^U EP , (16) 

where g a bcd = 9ac9bd — 9ad9bc is the de Witt supermetric and r] abc d is the 4-dimensional permu- 
tation tensor. It follows that 

E ab = C acbd u c u d , (17) 

Hob = \^acdPbe d U e , (18) 

where e abc = r] abcd u d is the alternating tensor on the observer's 3-dimensional rest space. The 
electric part of the Weyl tensor plays the role of the tidal tensor associated with the Newtonian 
gravitational potential, while H ab is essential for the propagation of gravitational radiation. 
Given that the Weyl tensor vanishes in FRW spacetimes, E ab and H ab provide a covariant and 
gauge invariant description of perturbations in the gravitational field. The electric and magnetic 
components of C a b c d also support the different polarisation states of propagating gravitational 
radiation and obey evolution equations remarkably similar to Maxwell's formulae pip. In the 
presence of a weak magnetic field, the linear evolution of E ab and H ab is determined by the 
system 

E ab = -QE ab + ±eiU - i/i(l + w)a ab + curing (19) 
H ab = -@H ab - cm\E ab + |curin a6 , (20) 

0"afe = -|©<7ab - E ab + \ll ab + D( a U 6 ) , (21) 

IU = -fGIU, (22) 

supplemented by the constrains 

B b E ab = ±D a/ u + \D a H 2 - \e abc H b cuT\H c , (23) 

B b H ab = n(l + w)u a , (24) 

D b IU = e abc H b cm\H c - \V> a H 2 , (25) 

D b a ab = §D a 6 + curkj a . (26) 

Finally, one should keep in mind that to linear order 

H ab = curla ab + D( a u 6 ) . (27) 

4.3 Isolating the tensor perturbations 

In the absence of the magnetic field one isolates the pure tensor perturbations of E ab and H ab , 
namely the gravitational waves, by demanding that 

u a = = D a/ u , (28) 

at all times. The above constraints set all the linear scalar and vector perturbations to zero, while 
ensuring that the remaining tensor fields are all transverse (i.e. T) a E ab = D a H ab = D b fj a ;, = 0). 
In the magnetic presence, however, one needs to impose two additional constrains 

B a H 2 = = e abc H b cmlH c . (29) 
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In other words, the spatial gradients of the magnetic energy density vanish and the field is also 
force free. Restrictions ( p9|) ensure that constraints (|2^) also hold in the presence of the magnetic 
field (see Eqs. (37), (41) in p6|]). Together, conditions ( |29| ) imply that D a p = 0, ii a = 

and D a G = 0. In particular, the 3-gradients of the pressure vanish as a result of the barotropic 
equation of state. The vanishing of the 4-acceleration comes from the linearised momentum- 
density conservation law (see Eq. (87) in fT5|). Finally, the linear propagation equation of 
gradients in the matter density guarantees that D a @ = to first order (see Eqs. (90), (91) 
in |15[). The same sets of constrains also guarantee the transversality of all the associated 
tensor fields, as Eqs. (25)-(|26|) immediately verify. 



5 Linear magnetized gravitational waves 
5.1 Evolution of the wave anisotropy 

Having isolated the tensor modes, we drop the acceleration terms from Eqs. ( |2l| ) and (|2j 
Moreover, the spatial flatness of the background means that the first order relation H a b = curios 
translates into cur\H a b = —D 2 a a b, where D 2 = D a D a is the projected Laplacian operator. It 



follows that we can eliminate the magnetic Weyl tensor from Eq. (19) in favour of the shear and 
therefore reduce the total number of the equations by one. Thus, the linear evolution of the 
magnetized gravity waves is monitored by 

&ab = — §®Cafc — E ab + \l\-ab , (30) 

Eab = -eE a b + ^eU ab -^{l + w)a a b-B 2 a ab , (31) 

ti ab = -fen o6 . (32) 

According to Eqs. ( |30|) and (|3l]), the field acts as a source of gravitational wave perturbations. 
Indeed, even when a a b = = E a b initially, & a b , E a b / because of the magnetic presence. This 
is not surprising at all, given that magnetic fields are natural sources of anisotropic stresses. As 
we shall see later, however, the specific form of these stresses (i.e. the tension properties of the 
magnetic field lines) means that the field presence can also suppress gravity wave distortions. 

Confining to large scales we can ignore the Laplacian term in the right-hand side of Eq. ( |3l| ) 
and reduce the above given system to 



&ab 



— §©0"afc — E a b + \ II a fe , (33) 

E ab = -eE a b + ±en ab -±fi(l + w)a ab , (34) 

ri o6 = -fen o6 . (35) 

5.2 Evolution of the wave energy 

Consider the magnitudes a 2 = a a ba ab /2 and E 2 = E a bE ab /2 of the shear and the electric 
Weyl tensors respectively. Once the pure tensor perturbations have been isolated, these are the 
only scalars one can invariantly construct from a a b and E a b- Moreover, a 2 and E 2 provide a 
direct measure of the wave's energy density and amplitude^ Their large-scale evolution comes 

4 The energy density of gravitational radiation is determined by the pure tensor (i.e. the traceless and transverse 
- TT) part H^J of the metric perturbation (e.g. see ]§)) 

Cgw = 2 ^ ' ( 36 ) 
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by contracting Eqs. (^3|) and (|34|) with a a b and E a b respectively. In particular we find 



(a 2 )- = -|6(7 2 -X- \H 2 Y,, (38) 
(E 2 )- = -2GE 2 - i/x(l + w)X - \QH 2 £ , (39) 

with X = E a i,a ab , S = a a }>rfrf and £ = E a },r] a r] b , while ?? a is the constant unitary vector in the 
direction of the magnetic force lines. The above system closes with the following propagation 
equations 

X = -\QX -2E 2 -\H 2 £ - fi{l + w)a 2 - \®H 2 ^ , (40) 
t = -§G£-£-±# 2 , (41) 
£ = -G£ -\^{l + w)^-\QH 2 , (42) 

where the magnetic energy density simply redshifts with the expansion 

(H 2 )- = -±GH 2 . (43) 



According to Eqs. (p0|), ( pl| ) (or equivalently (J33|) , (34)) the magnetic effects on a a b and 
E a b propagate via Il a b = —H/ a Hf,\. Since both a a b and E a b are also trace-free tensors, the 
magnetic effects on a 2 and E 2 propagate through the last two terms in Eqs. (|38|), (|39|), namely 
via the contractions £ = Oabrfrj 3 and £ = E a bi] a r] b . The latter describe the "squeezing" and 
the "stretching" of the space, along r] a , that is caused by the propagating gravitational wave. 
Crucially, r/ a is the direction the magnetic tension acts along. Hence, after the scalar and the 
vector modes have been switched off, the only remaining linear magnetic effect comes from the 
tension properties of the field lines. 

The scalars S and £ are related to spatial curvature perturbations via the Gauss-Codacci 
equation. Contracting Eq. ( JlOj ) twice along rj a and keeping up to first order terms we arrive at 

= £ - ieS , (44) 

where 5i is the total spatial curvature perturbation along the magnetic direction.^] Recall that 
the total contribution of the field to K is zero (see Sec. 3.2), which explains the absence of any 
magnetic terms in the above. Also, given that all scalar perturbations have been switched off 
and that the background geometry is flat, Eq. (fi4|) contains no matter or expansion terms either. 

where a is the dimensionless scale factor, a,/3 = 1,2,3 and the dash indicates conformal time derivatives. In 
a comoving frame the shear of the fluid flow is a 3-tensor (i.e. ctoo = = ctoq)- Also, the pure tensor part 
of the covariantly defined shear is related to the traceless ransverse component of the metric perturbation by 

<p = a ( H l0)', with <t = a ~ 3 ( H TT)' ( see H> EH)- 0n usin S these relations E q- © S ives 

PGW ~ °TT ' (37) 

with a 2 = 

5 I n M | the scalars E and E were related to spatial curvature distortions by twice contracting the PSTF part of 
Eq. (p. along r\ a . Here we exploit the fact that the magnetic contribution to 3-curvature perturbations (along 
the direction of the field lines) vanishes, and relate E and £ to the total spatial curvature perturbation K. This 
choice does not alter the essence of the calculation, but allows for a more transparent discussion on the subtle 
role of spatial curvature perturbations. 
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6 Magnetic effects on gravitational waves 



6.1 The radiation epoch 

When radiation dominates w = 1/3, \i = 3 /4i 2 and © = 3/2i (see Eqs. 
according to Eq. (Q), the magnetic energy density drops as 

2 



(15)). Then, 



(45) 



relative to a comoving observer with u a = 5®Uq. On using the above result, we find that in the 
same epoch the energy of the waves, as they propagate through the magnetized radiation fluid, 
is governed by the equations 

(a 2 )' = -2t~ 
(E 2 )' = -3f 

which form a closed system with the set 



E 



X 



±ai~ 2 £, 



\t- 2 x 



\ar ?J £ . 



(46) 
(47) 



X' = -\t~ x X -2E 2 -r 2 a 2 -\aT 2 £ - |at" 3 S , (48) 

£' = -£- \at~ 2 , (49) 

£' = -ft _1 £- |i~ 2 £- iat" 3 , (50) 

where a = H^tf) and the dash indicates derivatives with respect to coordinate time. Note that 
the zero suffix corresponds to the initial time to. The system (fPj|)-(|5"0|) accepts the late-time 
(i.e. when t 3> to) solutions 

So 



E 1 



ol + AE 2 tl 



2^0*0 — 2 
2 X X { 



So 



to 



1 tt2 
6 H 



#0*0 



to 



'Ml / C 1 S 1 jr-2 \ TT 

5 t - o- «-Hn #1 



*0 



(51) 
(52) 



for the magnitudes of and E a ^ respectively (see Appendix). Moreover, for radiation the twice 
contracted Gauss-Codacci equation (see (|44|)) gives 

^0 = &0 h , 

2 to 

which exactly coincides with the gravitational wave contribution to the parentheses in Eqs. (|5] 
(|52|). On using the above, solutions (f)l|), (|52]) transform into 



(53) 



a% = 5 [^o + 4^o*o " 2^oto] 



and 



i^o 
~ 2 t 



I (^0 - \Hq) #0*0 > 
t2\ tt2 



(jSt - ±H Q ) Hq 



(54) 



(55) 



respectively. The quantities in brackets describe the non-magnetized case. One can easily verify 



this by comparing our solutions to those of the magnetic- free studies [ 22 

Results (|J) and (|55|) show that the field leaves the evolution rate of a 2 and E 2 unaffected, 
but modifies their magnitudes. Therefore, the overall magnetic impact depends entirely on the 
initial conditions and it is sensitive to the initial value of 3ft. The latter describes perturbations 
in the spatial curvature, along the lines of the field, caused by the passing gravitational wave. 
In principle ?R.q can take either positive or negative values. 
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6.2 The dust epoch 



During dust domination w = and the zero order equations guarantee that = 2/t and 
[i = 4/3i 2 . Accordingly, the magnetic energy density evolves as 



h 2 = m 



8/3 



(56) 



relative to a comoving observer. The evolution of a and E is monitored by the system 



(E 2 y 



■§rV 



x 



-At^E 2 - \t~ 2 X - I3r ll/Z £ . 



■ft- 1 * - 2E 2 - §i~V 



supplemented by the set 
X' = 

s' = 

£' = 

with (3 = H^t^ 3 . Similarly to the radiation case before, the above has the following late-time 
solution 



(57) 
(58) 

(59) 
(60) 
(61) 



a 



E 



k i a o + 1^0*0 - §<*b*o] I j J 



_9_ 

25 



2/3 



_9_ 

50 



Sn 



2 S 

~ 3 t 



6^0 



Hn Hnt 



2.2 



fc 



2/3 



8/3 



') [ c 2^0 
50 1^0-3^ 



r 



(62) 
(63) 



In the dust era the twice contracted Gauss-Codacci equation (H3) gives 



^o — £o — | — 



2 ^0 
*0 



which again coincides with the wave contribution to the parentheses in Eqs. (|62|), (|6 
the above we may recast Eqs. (62) and (|63|) into 

(to\ 2 ^ / ^ \ 2 /3 



and 



f? 2 — 9 

^ - 25 



^ + 5i^ 



8/3 



S5 (^o - g-Ho) #o I — 



8/3 



(64) 
On using 

(65) 
(66) 



respectively, where again the brackets describe the non- magnetized case [22, 23]. As with ra- 
diation before, the overall magnetic effect depends entirely on the initial conditions. Note 
the magneto-curvature terms in Eqs. fl54|), ( |55"|) and (|65|), (|66|). Qualitatively, this magneto- 
geometrical effect tends to reduce the energy of the wave when the initial curvature distortion 
is positive (i.e. for 3Jq > 0), but increases both a 2 and E 2 if 3^o < 0. Quantitatively, the effect 
gets stronger with increasing curvature distortion. We will return to this non-trivial behavior of 
the field later. 
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6.3 Production of gravitational waves 

Consider the dust dominated era and assume that there are no gravitational waves originally 
present. In this case the initial conditions are a\ = E 2 = X§ = £$ = So = 0. Then solutions 
(H) and flU) give 

° 2 = T^(y) 2/3 , (67) 

E 2 = ^ 4 (^) 8/3 , (68) 

with an analogous result for dust. Thus, the magnetic presence has led to gravity wave pertur- 
bations. Note that this magnetically induced production of gravitational waves results purely 
from the anisotropic nature of the field. Qualitatively, the tension properties of the magnetic 
field lines are of no consequence. That is to say any other anisotropic source would have also 
triggered gravitational wave distortions. Indeed, results (|S7|), ( |68[ ) remain unchanged when the 
minus sign of the magnetic terms in Eqs. (I3§l)-([l2l) is replaced by a plus, namely when we switch 
from tension to ordinary positive pressure. Assuming the presence of a cosmological magnetic 
field at recombination we may use result (37) to estimate the strength of the induced grav- 



itational wave. Adopting an upper limit of 10~ 9 G, in today's values (see p9|1 ), we find an 
induced wave of the order of 10~ 70 GeV 4 also in today's values. The latter lies twelve orders of 
magnitude bellow the energy density of the inducing magnetic field. 

6.4 Boosting and damping of gravitational waves 

Let us now return to the magneto-curvature terms in Eqs. (fjl), (|55| ) and (|65|), (|66|). Their 
presence allows for the possibility of a zero overall magnetic effect if 3ft = H 2 /6 initially, namely 
if the curvature distortion along the field lines "equals" the isotropic magnetic pressure. Alone, 
the magneto-curvature terms reduce the energy of the wave if 5fto > 0, but lead to an increase 
when 3?o < 0. The effect results directly from the tension properties of the magnetic force lines. 



In fact, if tension were replaced by ordinary pressure in Eqs. (|3q)-(42) the effect is reversed. 
This intricate behaviour can be seen as the field's reaction to spatial curvature deformations. 
More specifically, the non-linear Gauss-Codacci equation ([!(]) gives 

K = n abV a V b = la 2 - ±ea abV a V b + Vc( a v C b)V a V b + , (69) 

for the total perturbation in the spatial curvature along the direction of the field lines. Note 
the absence of any magnetic terms in the above since the contribution of the field to 'R- a b r l a " r f ) 1S 
zero (see Sec. 3.2). Also, given that both matter and volume expansion perturbations have been 
switched off (see Sec. 3.3), 2/i — 2 /3 = due to the background flatness. The last three terms 
in the right-hand side of Eq. ( p9| ) describe the "stretching" and the "squeezing" of the space 
(in the direction of the field lines) caused by the passing gravitational wave and take positive 
or negative values. However on average, say over one oscillation period, their contribution to 
TZ a bi] a r] b amounts to zero, leaving a 2 as the only wave input to curvature deformations. Thus, by 
increasing a 2 when the initial curvature is negative and by decreasing it when 3f?o > 0, the field 
tends to minimise the curvature perturbation along its own direction. This magnetic reaction 
to curvature distortions is indicative of the tension properties of the field lines, that is of their 
tendency to remain as straight as possible. Analogous magneto-curvature effects have also been 
identified on the expansion of magnetised cosmologies (see JITJ] ) and might be interpreted as an 
indication of a magnetic preference for spatial flatness [pp| , [21[| . 
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7 Discussion 



The question of how magnetic fields interact with gravitational radiation is as timely as ever in 
view of the forthcoming gravity wave detection experiments and the ubiquity of magnetic fields 
in the universe. The lack of extensive research on the subject and the unique features of magnetic 
fields make the outcome of any such study difficult to foresee, while it could also probe unknown 
as yet aspects of fundamental physics. In the present article we have considered the impact of a 
large-scale magnetic field on gravity waves in the cosmological context. Our starting point was a 
spatially flat FRW background which was subsequently perturbed by weak gravitational waves 
and a weak magnetic field. We adopted a geometrical approach and employed the covariant 
formalism to examine the field effects on a 2 and E 2 , the scalars that directly describe the 
energy density of gravitational radiation. Throughout the analysis we maintained the pure 
tensor nature of the perturbed variables, employed scalars that were invariantly constructed 
from these tensors, and did not assume any a priori relation between the magnetic and the 
gravitational wave anisotropies. The geometrical nature of our approach brought to the fore the 
tension properties of magnetic fields and revealed their subtle interconnection with the spatial 
geometry of the magnetised spacetime. We found that the overall impact of the field depends on 
the initial set up and is particularly sensitive to the initial curvature deformation as measured 
along the direction of the field. In the absence of gravitational waves the magnetic presence led 
to wave production. Given the generically anisotropic nature of magnetic fields this is not a 
surprising result. The non-trivial effects came from the intricate coupling between the geometry 
and the tension properties of the field. The presence of the field was found to suppress the 
energy of gravitational waves when the initial curvature perturbation was positive, but led to 
a boost in the case of positive curvature deformation. Overall, the field seemed to react to the 
curvature distortion caused by the propagating wave and tried to keep it down to a minimum 
by modulating the wave's energy accordingly. 

The complete dependence of the magnetic effects on the initial conditions is rather unfortu- 
nate, as at this stage it is not clear what are the best physically motivated initial configurations. 
The ambiguity stems from the fact that, to linear order, 5?o depends on So and £q (see Eqs. (|53|) , 
fl64|)) which can take either positive or negative values. Ii is likely that a non-linear analysis 
would incorporate a 2 , namely the wave's original energy, thus favouring one initial configuration 
at the expense of the other. That aside, some interesting questions emerge when one is allowed 
to speculate on the basis of the linear results. Astrophysical magnetic fields, for example, are 
quite widespread and also considerably strong. Typical spiral galaxies carry extensive fields of 
few [iG and compact stars can locally support magnetic fields as strong as 10 16 G. If the field 
presence were to boost gravity wave perturbations in general, then, depending on the efficiency 
of the mechanism of course, one might think that detecting gravitational waves should have 
been a rather straightforward task. If, on the other hand, magnetic fields can suppress gravity 
waves, the ubiquity of cosmic magnetism could prove a considerable setback for the forthcoming 
gravity wave detection projects. 
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Appendix: Solutions 

One obtains the solutions (|5l|), fl52| ) and (|62|), (p3|) after a straightforward but rather tedious 
calculation. Here, we guide the interested reader through the intermediate steps. To begin with, 
Eqs. (|46|), ( |47| ) and the set (^)-(|50|) accept the solutions 



a 2 = _l Cl+C2t -3_ 2C3t -3/2_4 C4t -l + 4 C ^-5/2 + 4 a 2^ 2) (?()) 

£ 2 = -id^ + i^+Cat-^ + iC^ + iCgt^ + ^aV 4 , (71) 

and 

X = dt- 1 + C 2 t~ 4 + C 3 t" 5/2 + C A t~ 2 + C 5 t" 7/2 + |a 2 t~ 3 . (72) 

E = faT^ - fa-^st -372 - fa" 1 *" 1 , (73) 

£ = -l a ~ l C A r l - \aC^ /2 - \ar 2 , (74) 

respectively, where a = H^t^ and C t (with i = 1 , . . . ,5) are the integration constants. Solutions 



(|74|), (|73|) immediately provide the expressions 

C A = — ^ (2£ — E i 1 — f-^o) -^0*0 > (75) 
C 5 = - (5 + E t 1 + f Ho) H$t 9 /2 . (76) 



which substituted into Eqs. (72)-d74) lead, after a lengthy but straightforward calculation, to 



the rest of the integration constants 



C\ — — I [2Eq*q + \o\ - X t - (£0 - |E t 1 - \Hl) i?o*o] > (77) 
C2 = § [-#0*0 + °o + ^0*0 + (f £0 + §Eo^o 1 + if-^o) -^0*0] *o > (78) 
C 3 = ± [4£ 2 i 2 - 2<rg + Xoto + (4£o - l^oto 1 - ^H 2 ) #0*0] ijf . (79) 



On using results (|F5|)-(|79|), we arrive at the late-time (i.e. for t 3> to) solutions ( pip and ( [52] ) 
respectively. 

For the dust era we start from the system (|5"7| )-(|61~D and then proceed in a completely 
analogous way to obtain the late-time solutions (|65| ) and (|66|). Here, we only provide the full 
solutions to the set (|57|)-(|6"l|). They respectively are 

a 2 = 3 Clt -4_l C2t -2/3_3 C3t -7/3_ 2C4t -2 + 6 C5t -ll/3 + 3 ^ (g0) 

E 2 = ! Cl t- 6 -iC 2 t- 8 / 3 + 2C 3 t- 13 / 3 + ^ (81) 

X = Cit" 5 + C 2 t~ 5 / 3 + C 3 t- 10 / 3 + C 4 t~ 3 + Cgt" 14 / 3 + 3/3V 13 / 3 , (82) 

E = |^- 1 C 4 t- 1/3 - f /T^st" 2 - 2/3t" 5/3 , (83) 

£ = -l/r^t" 4 / 3 - ^/r^st" 3 - 2(3t- 8/3 , (84) 

where /3 = Hq^ 3 and C; (with i = 1 , . . . ,5) are the integration constants. 
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